Causal signal transmission by quantum fields. 
IV: The causal Wick theorem. 
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Wick's theorem in the Schwinger-Perel-Keldysh closed-time-loop formalism is written in a form 
where the place of contractions is taken by the linear response function of the field. This result 
demonstrates that the physical information supplied by Wick's theorem for operators is propagation 
of the free field in space and time. 
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I. INTRODUCTION 

In the previous papers of these series [lHl] we investi- 
gated response properties of quantised fields. Analyses in 
those papers were, so to speak, pure kinematics. Start- 
ing from this paper we extend our analyses to dynamics. 
Here we derive the formal tool on which other dynami- 
cal results depend, and which we term the causal Wick 
theorem. The latter is a functional (Hori's) [J, [5| form 
of Wick's theorem 6 8] in the Schwinger-Perel-Keldysh 
closed-time-loop formalism 9l4ll|. written in terms of the 
linear response function [12 1 characteristic of the field in 
question. The term causal Wick theorem was introduced 
in the unpublished paper [l3j . 

In the standard Feynman-Dyson technique of quantum 
field theory 0| , Wick's theorem supplies dynamical infor- 
mation about free fields, expressed by their contractions, 
or Feynman propagators (we avoid the term causal Green 
function, because we associate "causal" with "retarded," 
as in causal response and causal space-time evolution). 
In the closed-time-loop formalism, propagators becomes 
matrices, but the fact that Wick's theorem is a supplier 
of dynamical information contained in propagators does 
not change. 

As was shown in Refs. (see also 0,[H), 

compo- 
nents of the Perel-Keldysh matrix propagator may be ex- 
pressed by the linear response function (known in quan- 
tum field theory as the retarded Green function 0, 
characteristic of the field. As a physical quantity, this 
function emerges in Kubo's linear response theory (l2l | 
for the field in question, see Refs. 0, Q and sections ttV CI 
and III Dl of this paper. Linear response of the free elec- 
tromagnetic field was considered by Schwinger [16]. The 
only information about the field needed for Wick's the- 
orem is its response. Thus the formal operation of op- 
erator reordering which is the subject of Wick's theorem 
is inherently connected to how the field propagates in 
space-time. One may say that propagation of the field 
is the actual physical content of Wick's theorem. The 
causal Wick theorem makes this evident. 

An immediate word of caution is in place here. In the 
literature, the term Wick's theorem is used for two differ- 



ent facts. Wick's original theorem expresses time-ordered 
products of free-field operators as linear combinations of 
normally-ordered products of the same operators. This 
theorem holds irrespective of the quantum state of the 
system and ipso facto may be used for systems out of 
thermal equilibrium. The other meaning of this term is 
a Gaussian factorisation of Green functions of free fields 
in thermal equilibrium (underlying, e.g., the Matsubara 
diagram techniques 17]). Phase-space methods may be 
applied to quantum dynamics in both cases, cf. Refs. [lH 
and [ill [lj| , respectively. This particular paper is part 
of the effort directed at extending the out-of-equilibrium 
methods of Ref. [3] to space-time evolution of arbitrary 
interacting quantum fields, cf. the introduction to our pa- 
per [l| . The question of how far the equilibrium phase- 
space techniques of Refs. [H, Ql| may be generalised is 
outside the scope of this paper (and very much outside 
the scope of this series of papers). Consequently, here 
the term Wick 's theorem is used strictly in the sense of a 
relation between time-ordered and normally-ordered op- 
erator products. Had the need arisen to emphasise the 
said distinction, we shall talk about Wick 's theorem for 
operators. 

This paper is self-contained and does not directly rely 
on our previous work. A brief overview of papers [l|- 
Q seems nonetheless beneficial for putting this one in 
perspective. In Refs. [l|-[3| we demonstrated a structural 
equivalence of response properties of stochastic c-number 
and quantised fields. It applies without reservation to 
an arbitrary bosonic field, whether real or complex, free 
or interacting, in or out of equilibrium. At the cost of 
exploiting essentially nonclassical objects, such as Grass- 
mann sources [1, [H, [2(| , it extends also to fermionic fields. 
The technical tool making this equivalence evident is re- 
sponse transformation [Jul], cf. also (T3.ITH. Formally, it 
is a change of variables (termed response substitution) in 
the generating functionals of quantum averages of prod- 
ucts of Heisenberg field operators ordered in the closed- 
time-loop style. In [14j , a similar substitution was termed 
a change to causal variables. In this paper we use both 
terms: response substitutions introduce causal variables. 

In [lHU, causal variables were identified for the har- 



2 



monic oscillator, and then postulated for interacting 
fields. Papers [U, S| were concerned with consistency 
of such generalisation. It was shown that formal quan- 
tum response structures, which emerge as a result of the 
response transformation, coincide with typical response 
structures characteristic of stochastic c-number fields. In 
particular, classical reality and causality properties were 
shown to hold, exactly and without reservation, for quan- 
tised fields. 

The result of this and forthcoming papers in a nut- 
shell is that formal response structures identified in [H-Q , 
firstly, naturally extend to underlying quantum dynam- 
ics, and, secondly, under macroscopic (also called meso- 
scopic [28j |) conditions turn into physical response struc- 
tures. The causal Wick theorem allows one to rewrite 
Dyson's standard perturbative approach of quantum field 
theory directly in the response representation. Similarly 
to how response transformations reveal the kinematical 
response properties of quantised fields, the causal Wick 
theorem reveals the dynamical response properties hid- 
den in quantum equations of motion. For examples see 
Refs. PS [11112 • 

The causal Wick theorem may be seen as a step to- 
wards dispelling the myth about inconsistency of space- 
time evolution with relativistic quantum field theory. 
The reason why the Feynman-Dyson type of approaches 
do not allow one to describe propagation of quantum 
fields in space and time is purely technical: they do not 
operate with sufficient number of Green functions. This 
becomes obvious already if looking at pair Green func- 
tions. In the Feynman-Dyson techniques, there exists 
only one pair Green function per field — the dressed 
Feynman propagator. That this is not enough is obvious 
from a classical analogy. For a classical stochactic field, 
we have two pair quantities: the linear response func- 
tion and the average of a product of two fields. In the 
quantum response viewpoint, the corresponding quan- 
tities are Kubo's linear response function |12fl and the 
time- normal average of a pair of operators 0, HH, [23j | . 
One dressed Feynman propagator cannot be equivalent 
to two independent quantities. 

In the Schwinger-Perel-Keldysh closed-time-loop for- 
malism 0-[ll|) the number of independent pair Green 
functions doubles. Consequently, this formalism and the 
response representation are connected by a one-to-one 
transformation. This is all that response tranformations 
are about. Neither these transformations nor the causal 
Wick theorem may be introduced in the Feynman-Dyson 
techniques. Both are strictly subject to the closed-time- 
loop framework. 

We stress once again that we are talking about arbi- 
trary noncquilibrium systems where fluctiations and re- 
sponse are not generally related to each other. In vacuum 
the time-normal average vanishes, while in equilibrium it 
is related to response by Kubo's fluctuation-dissipation 
theorem. In both cases one enjoys the luxury of having 
only one independent pair Green function to work with. 
This underlies success of the standard Feynman-Dyson 



techniques, and of the so-called Retarded- Advanced for- 
malism [24, 25] , concerned with field propagation in inter- 
acting thermal systems 29]. Our results, we repeat, ap- 
ply to arbitrary nonequilibrium systems, where the said 
luxury is unavailable. 

We make extensive use of the functional (Hori's) form 
of Wick's theorem 0, @ (the Hori-Wick theorem, for 
short). Outside the quantum- field-theoretical commu- 
nity, Wick's theorem is often regarded as something high- 
minded, and its functional form as high-minded in the 
extreme. In fact, the Hori-Wick theorem is a major sim- 
plification. Firstly, it is an excellent mnemonic tool. Its 
algebraic rigidity leaves one only with the freedom of for- 
mulae for contractions. These are worked out applying 
Wick's theorem to pair products, which is next to trivial 
(cf. sections III Bl and II V Dl) . Secondly, it reduces compli- 
cated calculations to straighforward algebra (cf. sections 
III CI and HVEI) . Last but not least, the Hori-Wick theo- 
rem greatly simplifies extention of the formal techniques 
to fermions. By making use of auxiliary Grassmann vari- 
ables [IHHil, sig n rules characteristic of Wick's theorem 
for fermions are observed automatically, and the algebra 
remains mostly as for bosons. 

The paper is structured as follows. We start from de- 
riving the causal Wick theorem for the harmonic oscilla- 
tor in section |TTJ In section Hill we extend it to the real 
(Hermitian) field. In both cases, results of Ref. [2 form 
a coherent background for analyses in this paper. Un- 
like the real field, formulae we need for complex fields 
are absent in [1, 3]. Filling this gap is the subject of 
section IIV1 which is a sizeable fraction of the formal ef- 
fort in the paper. The causal Wick theorem for com- 
plex bosonic and fermionic fields is derived in section fVl 
In section I VII we outline simplifications characteristic of 
the nonrelativistic (resonance, rotating wave) approxima- 
tion. Appendices take care of formal particulars brushed 
under the carpet in the main body of the paper. They 
also outline details of derivations making the paper self- 
contained. Of independent interest to the reader may be 
appendices [B] and [C] which deal with the mathematical 
background. In appendix [B] we discuss variational deriva- 
tives by Grassmann fields. It is a somewhat streamlined 
version of the approach of Ref. Q . As a crucial generali- 
sation, it includes chain rules for variational Grassmann 
derivatives, which are shown to be identical to those for 
conventional variational derivatives. In appendix [Cl we 
show that Wick's theorem may always be written as a 
Hori-Wick theorem; in a sense, this statement is more 
general than Wick's theorem itself. 



II. THE CAUSAL WICK THEOREM FOR THE 
HARMONIC OSCILLATOR 

A. Hori's form of Wick's theorem 

To better orient the reader we start from a simple ex- 
ample. Consider an harmonic oscillator with the Hamil- 
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FIG. 1: The C-contour. 



tonian in the Schrodinger picture, 

Ha = fojJoa 'a, (1) 
where a^, a are the bosonic creation/annihilation pair, 

[a,at]=i. (2) 

In what follows we use condensed notation, 



f9 = / dtf(t)g(t), 



fKg 



dtdt'f(t)K(t - t')g(t'), 



(3) 



where K(t — t') is a c-number kernel, and f(t), g(t) may 
be c-numbers, q-numbers, or functional derivatives. Ex- 
cept for c-numbers, the order of factors in §3§ matters. 
Notation ([3]) saves space and makes formulae much more 
transparent. It also greatly simplifies extention to quan- 
tised fields, cf. section UTT1 below. 

We assume that the reader is familiar with the concepts 
of normal ordering, denoted and of the Schwinger- 

Perel-Keldysh's closed-time-loop, or C-contour, operator 
ordering 9Hll|. denoted Tc- In brief, normal ordering 
puts all aJ's to the left of all a's. The Tc ordering is best 
visualised as an ordering on the so-called C-contour de- 
picted in ng.[TJ Operators positioned on the forward and 
reverse branches of the C-contour are marked by pluses 
and minuses, respectively. The ordering rule coincides 
with the travelling rule indicated in fig. [TJ so that opera- 
tors on the forward (reverse) branches are time (reverse- 
time) ordered. For formal definitions see appendix I A II 

Now, let Q(t) be the time-dependent displacement op- 
erator, 



Qit) = ^-[a{t) + a\t)}, d{t) = 



-iujQt 



(4) 



Wick's theorem allows one to reorder Tc-ordered prod- 
ucts of many Q's normally. Its functional (Hori's [j, Q) 
form (the Hori-Wick theorem, for short) reads, 

TcHQ+^Q-) = -^n{Q):, where 

( 5 5 \ (5) 
T N (q + ) = exp Z C ^— , — j F{q + ,qJ)\ q _ =q+ . 

In these relations, J-(-, •) is an arbitrary c-number func- 
tional of two arguments, q±(t) are a pair of auxiliary 
c-number functional variables, and the reordering form 



Zc{-,-) is a functional bilinear form. In notation ([3]) it 
reads, 

= -j[v+Gfv+ - V-G* F V- + 2 V -G^ V+ ] , (6) 

where Gp and are the Perel-Keldysh contractions, 

(cf. endnote (30j ) 



- ihG F (t - t') = (0\T+Q{t)Q{t')\0) 

- ihG {+) (t - t') = (0\Q(t)Q(t')\0). 

The differentiations are performed using the formulae 
Sq+(t) 6q-(t) 



(7) 



6q+(t>) 5q-(t') 
Sq+(t) _ 6g-(t) 
5q-$) 5q+(t>) 



: S(t - t'), 

0. 



(8) 



Derivation of the Hori-Wick theorem for the oscillator 
may be found, e.g., in our Ref. appendix A. An alter- 
native derivation may be constructed using the method 
of appendix [I] which is readily adapted to a real field. 



B. "Mnemonic derivation" of the Hori-Wick 
theorem 

The Hori-Wick theorem ([5]) calls for few remarks, 
which should benefit the reader without a strong 
quantum-ficld-thcoretical background. As is shown in 
appendix [C] an exponentiated bilinear form of functional 
derivatives dutifully generates "all possible sets of con- 
tractions one can indicate, including the term without 
contractions" Q, irrespective of whether the "contrac- 
tions" make any physical sense. "Contractions" in this 
context are simply c-number kernels specifying the form. 
Hence mere awareness that Wick's theorem holds for a 
particular system allows one to immediatly write down 
the Hori-Wick theorem for it. It is convenient to firstly 
leave the reordering form unspecified, and then work it 
out applying the Hori-Wick theorem to pair products. 
For the oscillator, it means postulating eqs. ©,(16]), with 
Gf, G F and G( + ) regarded as a notation for three un- 
known kernels. They are found applying ([5]) to the prod- 
ucts, 



T c Q+(t)Q + (t') =T+Q(t)Q(t'), 
T c Q-(t)Q-(t')=T_Q(t)Q(t'), 
T c Q-(t)Q + (t') = Q(t)Q(t'). 



(9) 



The corresponding .F's are q+(t)q+(t'), q-(t)q-(t') and 
q_(t)q + (t'). Applied to pair products, the reordering ex- 
ponent in ([5]) may be reduced to two Taylor terms, 



cxpZc 



5q+ ' Sq- 



l + Z c 



5q+ ' 6q. 



(10) 
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FIG. 2: The three possible ways a pair of displacement oper- 
ators (bold dots) may be placed on the C-contour; Q — Q(t), 
Q' — Q(t'). Dashed lines signify the corresponding contrac- 
tions. 



It is easy to verify the formula, 



K(t-t'), 



(11) 



where K(t — t') is an arbitrary c-number kernel and 
f{t),g(t) are two auxiliary c-number functions. Using 
it we obtain, 

T+Q(t)Q(t') = -ihG F {t - t') + :Q(t)Q{t'):, 
T_Q[t)Q[t') = ihG* F {t' -t) + :Q{t)Q{t'):, (12) 
Q(t)Q(t') = -ihG^\t - t') + :Q(t)Q(t'):. 

Definitions (J7]) follow as a consequence of the relation, 

(0\:Q{t)Q{t'):\0) = 0. (13) 

Noteworthy is that three products ([§]) correspond to three 
possible ways of placing a pair of displacement operators 
on the C-contour (cf. fig. [5]), with in turn correspond 
to three unknown kernels specifying the form. Such 
"mnemonic derivations" allow one to quickly work out 
the details also in more complicated cases. 

C. The test-case formula 

The key formal results of [l| may be recovered by con- 
sidering the characteristic functional of vacuum averages 
of Tc-ordered products of Q(t), 

$v»c(»7+,»7-) = (0|T c exp (ir)+Q + -irj-Q-) |0>. (14) 

We use notation ([3]). This quantity expresses vacuum 
fluctuations of the oscillator. It may be explicitly calcu- 
lated using the Hori-Wick theorem ([5]), 



= expZ c - — , - — 
\oq+ oq- 



exp I 



</±=0 



ex.pZ c (ir)+,-ir)-), (15) 



cf. eq. (37) in |l| (heed should also be payed to endnotes 
31]). This relation as well as similar formulae below 
will be called test-case formulae. The intermediate ex- 
pression in (|15j) follows from ([5]), and from the fact that, 
for any J 7 ^, 



(0|:JV(Q) : |0> = F N (0). 



(16) 



A proof of the final one may be constructed using the 
relation [l5| . 

• Fl (^)- F2(5) l 9 =o=- F2 (^)- Fl(/) l/=o' (17) 

where f(t), g(t) are c-number functions and J-'i(f), ^(g) 
are c-number functionals. It can be verified expanding 
J~\ t 2 in functional Taylor series. Applying it to (|63p we 
find, 

$vac(ry+,?7-) = exp (iV+jjrJ exp ^ - ^-^r) 

xex P Zc{r 1 +,v'-)\ v , ±=0 , (18) 

where ri'±{t) are auxiliary c-number functions. The ex- 
ponents with derivatives here are nothing but functional 
shift operators. Applying the shifts and setting r]'±(t) to 
zero we recover eq. (|15[) . 



D. Response transformation revisited 

Response transformation [l[ relates contractions ([7]) to 
the retarded Green function of the oscillator, (cf. endnote 
301) 



G F (i - f) = G ( +\t - t') + G^\t' - t), 
G^(t - = G { + } (t - t') - G { k \t' - t), 



(19) 



where ^ stands for separation of the frequency-positive 
and negative parts of a function (cf. appendix I A 21) . The 
kernel Gr is Kubo's linear response function, 



G R (t - 1') = 



s(Q(t)) 

SJ c (t>) L=o 



= U(t-t')[Q(t),Q(t')], (20) 



defined with respect to the interaction with a c-number 
source J e (t), 



Hi(t) = -Q{t)J e {i) 



(21) 



with Q(t) being the Heisenberg displacement operator 
according to the Hamiltonian 



H(t)=H + Hj{t). 



(22) 



Explicitly, 



G R (t-t>) = e(t-t>) Skl " o{t - t ' ) . (23) 

We stress that we introduce one retarded transfer 
(Green) function, rather than a retarded and an advanced 
Green functions, as, e.g., in Ref. 24 1. 
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Equations (fl9)) suggest new variables r](x,t),j e (x,t) 
(causal variables, in terminology of 14]), such that, in 
notation (|3]), 

Zc(iV+, = ivGnjc- (24) 

The corresponding response substitution is [![, 



rj± = — ± 7^ +; , 



with the inverse substitution being, 

y = y+-y-, ie = + ?7-~ } 



(25) 



(26) 



By omitting the arguments we will be able to reuse these 
relations literally for a real bosonic field. We thus obtain 
vacuum fluctuations in an astonishingly classical form, 



$vac(»7+,»7-)| c . v . = eX P ( ir l G R3e), 



(27) 



where c.v. (short for causal variables) refers to substitu- 
tion (|25|) . 

With ([25j) being an integral transformation that smears 
rj(x,t) all over the time axis, restoration of classical 
causality in (|2T[) emerges as something of a miracle. 
Physical implications of a mere possibility of transfor- 
mation (|25p were discussed at length in Refs. PHI, [Hj] . 



E. Response transformation of the test-case 
formula 

So far our discussion stays within the results of Ref. pj . 
The step we undertake in this paper is to extend response 
transformations to Wick's theorem, which in turn allows 
one to extend them further to nonlinear quantum dy- 
namics. 

Formally, we want another pair of causal variables, 
C{t),Qe[t), to replace q±(t), such that eq. (fT5")) in new 
variables would read, 



4\ 



- exp 



3 5 

^G R — ) exp (irjq e + iO' ) | c=9o=0 



8q e 5( 



= exp (ir]G R je 



(28) 



We continue using notation ([3]). These variables may be 
identified comparing the linear forms in (1151) and in 



V+Q+ -V-Q- = VQe + Cje 



(29) 



Using (|23|) and the properties of the w operation (cf. eq. 
(|A10|) in appendix |A2|) we find, 



y+Q+ - y-q- 



= Je 



14 



■7° (+) 
J 11 



(30) 



This is satisfied by, 
'1+ 



q- 



q ± =q c ±K^. 



q { + +) 



,(-) 



(31) 



The new variables assure correct transformation of the 
functional differential form, 



5q+ ' 6q- 



dq e d( 



(32) 



Indeed, we know from Ref. [1] that eq. (pM)) is a con- 
sequence of transformation (jT9")) and substitution (|25l) . 
Now note that eqs. (pMf and (15"2"j) coincide up to the re- 
placements, 



»7± 



5 



Je 



5g±' JC ~5(' '' 
These replacements transform substitution (I25[) into, 

<5 id /^ (T) 



(33) 



i $ 



(34) 



Hence eq. (|32l) is found as a consequence of transforma- 
tion (01} and substitution (|3"3|). To show eq. ([22]) it suf- 
fices to verify eqs. (|34"|) . 

By the standard chain laws, 



Sq±(t) 



dt' 



saf) 



Sq e (t') 



5q ± (t)5((t>) 5q±(t) 6q e (t>) 



Using (|3"Tj) we obtain, 

sew i . ftzcfto (±) . 



(35) 



(36) 



where (^^(t) are the frequency-positive and negative 
parts of the delta-function (given by eq. (|A9[) in appendix 
IA 2[) . The second formula here is found writing the ex- 
pression for q D (t) explicitly as an integral transformation, 



9e(*) 



dt'[6^(t-t')q + (t') 

+ S(-\t-t')q_(t')}, (37) 



cf. eqs. (|A8|) . ([A9]) in appendix [A2] Combining ((35]) and 
([55| and again recalling eqs. (|A8[) . (|A9|) we have, 



±- 



M*0 



1 

= ±- 



h5((t) 



(38) 



This coincides with (|34l) up to overall phase factors. 

Few words regarding the choice of factors in eqs. 
(I3"2l . To start with, Planck's constant is absent from the 
RHSs of eqs. ^ and (|2"8"]h The way fi is distributed 
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among the variables in (|3Tj) warrants that q e has the 
dimension of the "field" Q. Similarly, j e (x,t) has the 
dimension of field source. The explicit i factor in (|27|) 
assures consistency with conventional definitions of char- 
acteristic functions and functionals in probability theory 
and theory of random processes. When comparing re- 
sults of this pap er to Ref. [l|, heed should be payed to 
endnotes I30L kill. 



F. The causal Wick theorem for operators 

Transformation (1311) was derived for vacuum averages, 
but nothing prevents one from applying it directly to 
eq. ([5]). As may be seen from eq. (|3Tj) . the condition 
q+(t) = q-(t) in terms of the causal variables becomes, 



q+(t) = q_(t) = q e (t), ({t) = 0. 



(39) 



With this observation, and accounting for (f3"2"j) . we find 
the second of eqs. © in the form, 



J^N(q c 



Sq c 8( 



rxp| - , — )j 7 (q c + hC i ) ,9c-^C (+) )| c=0 - 



(40) 



The first of eqs. ((5|) remains unchanged. Unlike f|28j> . 
eq. (|4l7| holds irrespective of the quantum state of the 
system. 



the "mode sum" ^ K and of the "spatial integration" 
/ dx, is problem-specific. Nothing precludes k from being 
(partly) continuous, and x from being (partly) discrete. 

The case of a free bosonic field was amply covered in 
our Ref. Transfer function of the real field is given 
by Kubo's formula, 

G n (x,x',t-t') = l -6{t-t')([Q{x,t),Q(x' ,t')]). (43) 

The Perel-Keldysh contractions read, 

-ihG F (x,x',t-t') = {0\T+Q(x,t)Q{x',t')\0), 



-ihG (+) (x,x',t-t') = (0\Q{x,t)Q(x',t')\0). 



(44) 



Explicit expressions for the kernels are found using (|4ip . 
Relations for the contractions in terms of the transfer 
function were worked out in Ref. QL 



G F (x, x', t-t') = G^ +) (x, x', t - t') 

-f" Cr (x , X) t £), 

G* + '(i, x', t - t') = G { + } (x, x', t - t') 

-G { ~\x',x,t' -t). 



(45) 



Equations (|43| . (j44)) and (|45|) generalise, respectively, 
eqs. © and tfTHj) . 



B. The causal Wick theorem for the real field 



Equation (|4T)f and its couterparts derived below are 
the main result of this paper. In what follows, eq. (|40)) 
will be extended to quantum fields in the true meaning 
of the word including fermions. 



III. EXTENSION TO THE REAL (HERMITIAN) 
FIELD 

A. Definitions and recollections 

Since response transformations apply mode-wise, the 
case of the real field is a trivial generalisation of that of 
the harmonic oscillator. The free-field operator is, 



Q(x, t)=Y, \l ^7 M*K e ~^' + <{x)die^ 1 ] , 



(41) 



where a£.,d K are bosonic creation-annihilation pairs, 

[a K ,dl]=S KK ,. (42) 

The mode frequencies uj k and the c-number mode func- 
tions u K (x) may be arbitrary. The meaning of the "co- 
ordinate" x and of the "mode index" well as of 



Similarly to the case of the oscillator we make use of 
the condensed notation, 



(46) 



fg = / dxdtf(x,t)g(x,t), 



fKg — / dxdx dtdt f(x, t)K(x, x ,t — t )g{x , t ) 



where K(x,x',t — t') is a c-number kernel, and f(x,t), 
g(x,t) may be q-numbers, c-numbers or variational 
derivatives. This notation applies throughout the rest of 
the paper. For the purposes of section [IV] below, f(x,t) 
and g(x, t) may also signify Grassmann fields and Grass- 
mann derivatives (see section IIV C 21 and appendix |5J . 
Except for c-numbers, the order of factors in (|46p mat- 
ters. 

Understood in line with (|4"6")l . all critical relations for 
the harmonic oscillator apply literally to the Hcrmitian 
field. This implies that eqs. (|2"3"]). © and $F5§ are 
substituted by l[4"5|), (|44]) and (gSJ), and that all func- 
tional variables acquire the second (spatial) argument, 
?7± (£) — > r)±(x,t), q±(t) —> q±(x,t), etc. Since the W 
operations apply only to the time variable, the spatial 
arguments are carried through all manipulations without 
affecting them. In particular, this applies to the trans- 
formation of the linear form (|29[) on which the derivation 
hinges. Response substitutions (|25p and (l3~Tj) . which we 
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have foresightedly written without arguments, apply lit- 
erally. 

Here we briefly reiterate the logic of the argument 
adapted to the real-field case. 

1. Formulate the Hori-Wick theorem, eq. ([5]). This 
implies the definition of the reordering form by eq. 
([6]), subject to notation (|46|) . and of the contrac- 
tions by eq. flUJ|. 

2. Use the Hori-Wick theorem to derive the test-case 
formula (fT5j) . 

3. Apply response transformations to the final expres- 
sion in the test-case formula, cf. eq. (|27|) . 

4. Work out the change of variables in the reordering 
form postulating transformation of the linear form 
(|29p . Verify transformation of the reordering form 
([521) b y proving eqs. (EH)) . 

5. Extend the causal Wick theorem to operators re- 
sulting in eq. (l4TI|) . 

Steps 1-3 of this schedule are a result of Refs. [H, H|, while 
steps 4 and 5 are the subject of this paper. In what 
follows this logic is further extended to complex fields, 
bosonic as well as fermionic. 



IV. TRANSFER PROPERTIES, WICK'S 
THEOREM AND VACUUM AVERAGES FOR A 
LINEAR QUANTUM CHANNEL 

A. Preliminary remarks 

In the rest of the paper, we extend analyses of sections 
HIl and Hill to complex fields, bosonic and fermionic. This 
takes a few generalisations of the results of Refs. [l], 
which are the subject of this section. These generalisa- 
tions may be divided into compulsory and elective. Com- 
pulsory generalisations are what in the introduction was 
called "filling the gap:" derivation of an analog of the 
test-case formula (fT5)) for complex fields, absent in [2,[3j]. 
The elective generalisation is an extension of complex 
fields to linear quantum channels. 

Quantum channel is a convenient formal way of mak- 
ing our analyses applicable to relativistic quantum fields 
with notrivial matrix structure. So, a conjugate pair of 
free fields in papers [l|, H[ was understood as a Hermitian- 
conjugate pair, whereas the Dirac-adjoint of the four- 
component spin-half field (say) involves change of sign 
of two components. To cover such cases we make co- 
efficient functions of the "straight" and "adjoint" fields 
independent, calling such pair a linear quantum chan- 
nel. Conventional quantised fields, whether bosonic or 
fermionic, are special cases of quantum channels. 

The algebra of Green functions of a quantum channel 
is in essence that of a complex field (see appendix [Dj . 
Results of Refs. 0, Q may therefore be used for quantum 



channels with minimal amendments. It is "filling the 
gap" that takes most of the formal effort in this section. 
Extension to quantum channels reduces to a large extent 
to the observation that "it does not change the algebra." 

B. The channel operators 

A linear quantum channel (linear channel, quantum 
channel, or just channel) is specified by a pair of "tilde- 
conjugates," 

(47) 

where u K (x),v K (x),v K (x),u K (x) are four independent c- 
number functions. The meaning of the "coordinate" x 
and of the "mode index" well as of the "mode sum" 
and of the "spatial integration" J dx, is problem- 
specific, cf. the remarks after eq. (|4T|) . The opera- 
tors b K , b' K ,c K , c' K form two standard creation/annihilation 
pairs for particles and antiparticles, which may be 
bosonic or fermionic, 

[Ml] T = [c«,4] T = <W. (48) 

Here, [•,•]— stands for a commutator and [•,•]+ for an 
anticommutator, 

[kX'} T = b K bl,Tbt,b K , (49) 

etc. In eqs. (|48|) . (j49]l. the upper sign is for bosons and 
the lower for fermions. To avoid lengthy specifications, 
we shall also use blanket terms commutator, commute, 
etc., with respect to both commutational and anticom- 
mutational relations. So, we say that, except (|48J1 . the 
creation and annihilation operators pairwise commute, 
meaning that they commute for bosons and anticommutc 
for fermions. 

The quantum channel is a natural generalisation of 
a complex quantised field, bosonic or fermionic. Equa- 
tion (jf?} obviously alludes to Dirac's fermions (zj. In 
this case, the tilde stands for the Dirac adjoint, x 
comprises the coordinate and the spinor index, k is a 
combination of the momentum and spinor index, and 
u K (x),v K (x),v K (x),u K (x) are suitable solutions of the 
Dirac equation. In general, these symbols may be as- 
signed any meaning. 

Definitions (|47[) are semirelativistic: they assume exis- 
tence of particles and antiparticles but not the relativistic 
covariance of the theory. The semirelativistic case in- 
cludes relativistic theories but does not reduce to them. 
Nonrelativistic solid-state theories dealing with particles 
and holes are in our sense semirelativistic. 
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C. Response problem for the quantum channel 

1. The Hamiltonian 



The pair of "modes" (|47|) becomes a channel when one 
couples to them a pair of external sources. Formally, this 
is expressed by the Hamiltonian, 

H(t) = H Q + H I (t), 

K 

Hi(t) = - I dxHi(x,t) 

dx[?p(x,t)a(x,t) + a(x,t)ip{x,t)\ . (50) 

By definition, this Hamiltonian is written in the inter- 
action picture, with Hq regarded free Hamiltonian, and 
Hi(t) interaction. The channel operators are interaction- 
picture operators. The corresponding Heisenberg opera- 
tors are &(x,t) and ^>(x,t). 

2. Summary of Grassmann fields 



For bosons, the sources in (|50|) are c-numbers, for 
fcrmions, they are Grassmann fields [3[ . Unlike Beresin 
and Vasil' ev [llll, we define Grassmann fields as linear 
combinations of Grassmann generators, and the deriva- 
tives by them as anticommuting variational derivatives. 
To make the latter unique, we regard dimensionality of 
the underlying Grassmann algebra a variable parameter 
For mathematical details see appendix [Bj 

For all practical purposes, Grassmann fields are formal 
quantities which commute with c-numbers and bosonic 
operators, and anticommute between themselves and 
with fermionic operators, e.g., 



cr(x, t)a(x' , t') — ±a(x' ,t')a(x,t), 
cr(x, t)a(x , t') = ±a(x' ,t')a(x,t), 

a(x,t)^(x' ,t') = ±-tp(x',t')a(x,t), 



(51) 



etc. As in (|48j) . the upper sign here is for bosons (when 
a and a are c-number fields), and the lower for fermions 
(when a and a are Grassmann fields) . Differentiation by 
Grassmann variables follows simple rules, 



SC 



= o, 



5a(x, t) 
S[X + y] 6X 



sy 



Sa(x,t) 8a{x,t) Sa(x,t) 
5[a(x',t')X] 



(52) 



Sa(x, t) 



= 5{t-t')8{x,x')X 



±a(x',t') 



SX 



where C is an arbitrary quantity which does not depend 
on a(x,t), and X, y are arbitrary quantities which may 
or may not depend on a(x,t). The last of eqs. (|52|) 
holds for c-numbers with plus and for Grassmann fields 
with minus. Furthermore, for linear substitutions with c- 
number coefficients, chain laws are identical for c-numbcr 
and Grassmann variational derivatives, cf. eq. (IB19[) in 
appendix lB 41 Together with the chain laws, eqs. ([5"Tj) and 
(1521) suffice for all algebraic manipulations in the paper. 
The rule of thumb is, watch the order of factors and how 
you commute/anticommute them, otherwise proceed as 
with c-numbers. 



3. Transfer functions of the channel 

We are interested in the linear response functions of 
the quantum channel, 



A R (x,x',t-t') 



A R (x',x,t' -t) = 



5<S/(x,t) 
Sa{x' ,t') I 'U=5-=o' 



6^(x,t) 



5a(x',t>) 



(53) 



They describe propagation of quantum signals in the 

channel from ip to tp and from ij) to tp, respectively. 
For a linear channel, these signals do not interfere with 
each other. That the derivative must be kept inside 
the average is a fermionic/Grassmann feature: Grass- 
mann derivatives do not commute with quantum averag- 
ing (Ref. [3, section 4.2). 

Similar to the harmonic oscillator, we introduce only 
retarded Green functions, rather than more commonly 
used retarded and advanced functions [24[ . It should not 
be overlooked that Ar is defined with a transposition of 
arguments. Such definition conveniently absorbs trans- 
positions which otherwise occur in response transforma- 
tions [H-Q. It also leads to natural algebraic structures 
in the case of fermions. 

Quantum dynamics with Grassmann sources was con- 
sidered, e.g., in Ref. H, section 5. What matters is that, 
so far as we preserve the order of factors in Hi, pertur- 
bation theory with this Hamiltonian may be constructed 
the usual way. It is then straightforward to obtain the 
explicit formulae, 



A R (x,x',t-t') 



-e(t-t')($(x,t),i>(x',t')]J, 



(54) 



and 



Sa(x, t) 



A K (x', x, t'-t) = --6(t - t'){ $(x', t'), ^{x, t)] T ), 

(55) 

cf. Ref. [3J . Both in (|54l) and (|5"5"|) the commutators are c- 
numbers, so that the quantum averaging is in fact redun- 
dant. The anomalous response functions of the channel 
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FIG. 3: The ten possible ways a pair of channel operators 
(bold dots) may be placed on the C-contour; j/i = ip(x,t), 
ip' — tp(x',t'), and similarly for the tilded operators. Dashed 
lines signify the corresponding contractions. Only those on 
the top row are nonzero. 



vanish, 



si?(x,t) 

8a(x',t') 



\a=a=0 



5^{x,t) 
Sa(x',t') 



cr=cf=0 



0. (56) 



Details of the calculation, which may be interesting as an 
example of algebraic manipulation of Grassmann fields, 
are consigned to appendix [El 

Up to transposition of arguments and the ±i/H fac- 
tors, the transfer functions are the retarded and advanced 
parts of the same field commutator, and must therefore 
be related to each other. Indeed, in a real problem, the 
said commutator obeys some sort of free equation. Know- 
ing its solution in one half of space-time suffices to find it 
in the other half. However, there does not seem to exist 
a closed relation between the two, valid for an arbitrary 
quantum channel. We treat Ar and Ar as if they were 
independent quantities. 



D. The Hori-Wick theorem for the quantum 
channel 



Definitions of the Tq and normal orderings for bosonic 
and fermionic operators were discussed at length in 
Ref. |3j, section 4.3; they are briefly reiterated in ap- 
pendix IA II The Hori-Wick theorem for the quantum 
channel reads, 



TcF{^ + ,ili + ,^-,i>_) = :J r N{'4 ) , ^)-, where 
( 8 8 8 8 



T N (ip + ,ip+) = expZ c 



\8ip+ ' S^p + ' St/)- ' (5^_ 



x^ + > + >_,^)|^ +jW+ . (57) 



In this formula, •, •, •) is an arbitrary c-number func- 
tional of four arguments. The reordering form in notation 



06) is, 

Zc(f+J+J-,f-) = -ie f h[f + A F f + - /_A F /_ 

+ /_aW/ + -/ + aH/_], (58) 

where 

-iHA F (x,x',t-t') = (0\T + ^(x,t)tjj(x',t')\0) 7 
ihA F (x,x',t-t') = (0\T^Tp(x,t)>4>(x' ,t')\0), 
-ihA {+) (x,x',t-t') = (0\<ip(x,t)i>(x',t')\0), 

ihA { - ] (x, x', t-t') = e f (0\i>(x', t')^{x, i)|fj), 

are the four contractions characteristic of the channel. 
We introduced the fermionic sign factor, 



(59) 



-1 for bosons, 
-1 for fermions. 



(60) 



The ± symbol is hereafter reserved for purposes other 
than distinguishing bosons from fermions, in particular, 
for distinguishing the C-contour indices as in eq. (p9"[) 
below; this notation applies throughout the rest of the 
paper. The ip's in (I5T1) and the /'s in ([551) are auxiliary 
functional variables, c-numbers for bosons or Grassmann 
fields for fermions. Algebraic rules for Grassmann fields 
were summarised in section TlVC 21 for their justification 
see appendix [Bl Proof of eq. (|57[) is given in appendix [Fl 
Similar to the oscillator, mere awareness that proof of 
the Hori-Wick theorem ([571) exists allows one to construct 
its "mnemonic derivation" in detailed analogy with sec- 
tion [IIBJ An arbitrary quadratic form of four functional 
arguments is specified by ten kernels, which correspond 
to ten possible ways of placing a pair of channel opera- 
tors on the C-contour, cf. fig. [3] Applying (|57|) to such 
pairs establishes the connection between kernels and con- 
tractions. This calculation utilizes a generalisation of eq. 
([TT]) to complex fields, 



±K-L)/U.I)fu'.l') 



e f K(x,x',t-t'), (61) 



where K(x,x,t — t') is a c-number kernel and 
f(x,t),f(x',t') are c-number functions (for bosons) or 
Grassmann fields (for fermions). Fermionic sign factors 
in ([55)1 and ([6"Tj) cancel each other. Of ten contractions, 
only the four on the top row in fig. [3] arc nonzero. They 
are given by ([59]) . 



E. The test-case formula 

As a test case when deriving response transformation 
of the Hori-Wick theorem we shall use a formula for the 
vacuum average, 

*vac(^+,?7+,^-,?7-) 



(0|Tc exp («?7 + -0 + + i'0+ 7 7+ — ifj-i>— — m/*_J7— ) |0)- 



(62) 
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As for the oscillator, the test-case formula expresses cal- 
culation of $ vac by means of the Hori-Wick theorem, 



expZ c 



5ip + Sip- Sip- 
x exp (irj + ip+ + i4>+v+ - ifj-ip- - | - =Q 

= exp [efZ c (ifj+ , in + , -if}-, -in-)]- (63) 

For bosons, the 77's and tp's are auxiliary c-number func- 
tions, for fermions, they are auxiliary Grassmann fields. 
The intermediate relation in (|63l) is found applying the 
Hori-Wick theorem to eq. (|62|) and using that, for any 
functional J-jv, 



(0|:J"jv(^):|0) - ?n(0,0). 



(64) 



The final formula is now subject to verification. 

For bosons, proof of (l63l) is a straightforward gener- 
alisation of that of eq. (fT8|) in section III CI The only 
difference is that we now have four shift operators to ap- 
ply. To extend (|6"5|) to Grassmann fields, we note that 
the sign factor in the exponent there cancels the "inbuilt" 
sign factor in Z Cl so that the last expression in eq. 
actually reads, 

®v ac (v+,V+,V-,V-) = exp{ih[fj + A F r] + - fj-A F r)- 

-^_A (+ V + r) + A ( ~V]}- (65) 

This formula has been verified for c-numbers. To extend 
it to Grassmann fields, it suffices to prove that, in both 
cases, 

$vac(*)+,f?+,'7-,»?-) 

= F(77 + A F 7 /+ ,7kA F 77_,77_A( + V,77 + A(-V), (66) 

where F(-, -,-,•) is a unique c-number function of four 
arguments. Equation (|66[) is proved in appendix [Gl We 
therefore regard the test-case formula (|63l) proved for 
both bosons and fermions. 



V. THE CAUSAL WICK THEOREM IN THE 
SEMIRELATIVISTIC CASE 



A. Response transformation of contractions and 
the causal variables 



As is shown in appendix Q2 the algebra of Green func- 
tions of a quantum channel is in essence that of a complex 
field, bosonic of fermionic (T), Q . We may therefore use 
results of those papers with minimal generalisations. So, 



(67) 



response transformations of contractions read, 

A i+ \x,x',t-t') = A R + \x,x',t-t') 

- A^ (x,x',t- t'), 
A^{x,x',t-t') = A < R ~ ) (x,x',t-t') 

-A ( -\x,x',t-t'), 
A F (x,x',t-t') = A R +) (x,x',t-t') 

+ A { R \x,x',t-t'), 
A P (x,x',t-t') = A { R \x,x',t - t') 

+ A R + \x,x',t-t'). 



These relations coincide with those derived in Refs. [l|, 
except for two minor differences: there are now two (for- 
mally) independent response functions, Ar and Ar, and 
Ar was defined with a transposition of arguments. 

Equations (|67|) hint at new functional variables, rj(x,t), 
fj(x,t), a e {x,t) and a e (x,t), such that, 

Z c (irj + ,ir) + , -ifj-, -if/-) = i£ f (r)A R cr c - ct c A r ?7). 

(68) 

We continue using condensed notation (|46| . The bilinear 
forms in (|68j) clearly allude to the response properties of 
the quantum channel, see section HV CI It should not be 
forgotten that evolution in Ar is from left to write, see 
eq. (|53p and ensuing comments. The causal variables are 
introduced by the response substitution [H-0], 



V± = ^TV iT) , V± = ^±V^\ (69) 
while the inverse substitution reads, 



fj = fj + - ?7_, (7 C = h[fj^' + fji_ } 



(70) 



To verify (|68p . it suffices to apply eqs. (|59j) and use the 
relations, 



fK^f = fKf& = f^Kf. 



(71) 



The latter follow from eqs. (|A10[) in appendix IA 21 For 
the reasons explained at the end of section III D\ we moved 
the imaginary unit from the transformations (|69|) , (f70|) to 
the quadratic forms . Otherwise our equations follow 
phase conventions of Ref . , cf . also endnote [3l[ . 



B. Response transformation of the test-case 
formula 



According to the schedule of section MI Bl we have 
completed the preliminary tasks 1-3, and are now ap- 
proaching task 4. We wish to find functional variables 
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t), ((x, t),ijj e (x, t), ip e (x, t) bringing the test-case for- 
mula to the form, 

^ va , c {fj+,V+, V-,V-)\ C , V , 

( ■ 6 A 6 \ ( ■ 6 X 5 

= exp - let— A R -~ exp i£ f — A R7T 

V He KJ V H He 

x exp (ifji> e + i(<r e ) exp ( - itp c i] - ia c () | - =f= £ =Q 
= exp [if)A R a c - ia c A R 7j] , (72) 



where c.v. refers to substitution (|69|) . We split exponents 
so as to emphasise that eq. (1721 in fact factorises. Fac- 
torisation is already seen in eqs. (j69l) . which express two 
nearly identical transformations of the "tilded" and "un- 
tilded" variables. Transformation of the former coincides 
with (f2"5)) , ^ofy , while that of the latter differs in the sign 
of 77. We therefore can recycle the algebra of section Hi El 
with minimal amendments. Postulating equality of the 
linear forms in (|63p and in (|72p we find the substitutions, 



*±=iTfi( (T) , ^ ± =^ c ±hC ( - T) , (73) 



where 



(74) 



is identical to (1311) , while that of the "untilded" ones dif- 
fers in the sign of Transformations laws for the deriva- 
tives by "tilded" variables are therefore identical to (j3"4"|) , 
while for the "untilded" ones we have to adapt (|3"4")l to 
the said sign change, 



Hi 



S 

H± 



i_5_ 

hSQ \Hc 
i S 



(T) 



(T) 



(75) 



In continuing analogy with the real field, these transfor- 
mations emulate substitution with the replacements, 



5 . <5 

Sip± H± 



5 

—1—, 

K 

5 

He 



(76) 



He 



so that the transformation of the quadratic forms of 
derivatives reduces in fact to eq. (|6"5|) . This concludes 
step 4 of the schedule of section IlIIB I for quantum chan- 
nels. 

C. The causal Wick theorem 



Similar to the real field, these transformations are 
rescaled versions of (|69l) . (f70|). Comparing them to eqs. 
(151) we see that transformation of the "tilded" variables 



Step 5 of the aforementioned schedule is straightfor- 
ward. Applying transformations (|T3[) to the formula for 
Tn in the Hori-Wick theorem (|57|) we find, 



exp 



6 S S - 6 

He S( H S^ e 



- ^e + fiC (_) , ^ + ^C (+) , A - hC {+) ) 1 c= . =o . (77) 
I 



The first of eqs. (|57|) remains unchanged. Equation ([77]) 
generalises eq. (1721) to operators and ipso facto to arbi- 
trary quantum states of the field (channel) . 



VI. THE CAUSAL WICK THEOREM IN THE 
NONRELATIVISTIC CASE 

A. Motivation 

The nonrclativistic case of quantum channel follows by 
dropping antiparticles (holes) in the channel operators 
(|47|) . that is, formally, by setting, 



v K (x) = V K (x) = 0. 



(78) 



This also covers a physically important case of quantum 
electrodynamics under the resonance, or rotating wave, 
approximation (RWA), commonly used in spectroscopy 
and quantum optics. Under the RWA, one works not 
with the real field operator (PUT) , but with its frequency- 
positive and negative parts, 



)(-)/ 



(79) 



K 

Up to notation, 

tp{x,t) = Q (+) (x,t), i(x,t) = Q { ~\x,t), 



(80) 
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eqs. (|79[) define a nonrelativistic quantum channel 
teractions (f2"Tj) generalised to the real field and (HT 
the channel coincide under the RWA, 



In- 
fer 



J dxQ(x,t)J e (x,t) 

dx [a(x, t)tp(x,t)+a(x,t)i>(x,t)], (81) 



subject to another instance of notation, 

a(x,t) = J<; +) (x,t), a(x,t) = Jt\x,t). (82) 

This leads to redefinition of the response function by the 
first of eqs. (J53J), and of the contractions by eqs. 

Heed should be taken of the fact that redefinition 
of the response function under the RWA implies that 
we do not regard a(x,t) and a(x,t) as being, respec- 
tively, frequency-positive and frequency-negative. In- 
deed, Kubo's theory assumes that no conditions are im- 
posed on the sources. However, with o~{x, t) and tr(x, t) 
constrained to the subspaces of frequency-positive and 
negative functions, we must also use differentiation with 
constrains. The latter is defined for an arbitrary func- 
tional F(-) as [HI, 

ST(Jr) 



Sjt\x,t) 



= / dt'S^(t'-t) 



ST J 



6F(J ) 
5J e (x,t') \_SJ e (x,t) 



l(T) 



(83) 



Irrespective of any particulars, such constrained differen- 
tiation cannot lead to causal (retarded) response because 
of the final separation of the frequency-positive and neg- 
ative parts. We return to the discussion elsewhere. 

A somewhat confusing feature of the nonrelativistic 
case is that, physically , it emerges as an approximation 
to more complicated cases, while, formally, it comprises 
a coherent set of exact relations. As was just illustrated, 
this exactness is a consequence of a consistent set of ap- 
proximations (see also a remark at the end of appendix 
B.4 in (3). A physically meaningful discussion of this 
point is postponed till particular models are approached. 
In this paper we treat the nonrelativistic case as a mere 
formal consequence of eq. ([751) . 



B. Response transformations of contractions for a 
nonrelativistic quantum channel 

Owing to (I75|) . the frequency-negative part of the com- 
mutator vanishes, 



A ( -\x,x',t-t') =0. 



The algebra of Green functions, expressed by eqs. (|D5[) 
(|D8|) . (|D10[) in appendix [Dl then becomes, 



A F (x, x', t-t') = A R (x, x', t - t') 

= 6{t-t')A {+ \x,x',t-t'), 

A F (x, x', t-t') = A R (x, x', t - t') 

= -9(t' -t)A {+ \x,x' ,t-t'), 

A i+) (x,x',t- t') = A R (x,x\t- t') - A R (a;,a; / ,t- t'). 

(85) 



Equations (|85| already contain the expression for Ap, Ap 
and A( + ' by Ar, and Ar. No additional transformation 
is necessary. 



C. The test-case formula. . . 

The nonrelativistic case is not covered in our 
Refs. p], |j| (it was considered in Ref. [3], but none of 
the formulae we need were derived there). The good 
news is that, formal differences notwithstanding, the key 
relations coincide in the semirelativistic and nonrelativis- 
tic case s. We therefore recycle as much as possible of 
the semirelativistic case, mostly marking the differences. 
So, the Hori-Wick theorem (|57|) and the test-case formu- 
lae ((55]), ([55)1 equally apply in the nonrelativistic case. 
We can also recycle eq. ([55)1 postulating causal variables. 
The definition of the causal variables is already different. 

Satisfying eq. (|68|) in the nonrelativistic case is 
straightforward. Insertion of eqs. (|85j)) into the LHS of 
(1581 yields, 

Z c (ifi + , irj + , -ifj-,-ifj-) 

= iefh[(f) + -ij-)A R r)+ -?7_A R (77_ -r?+)]. (86) 

Comparing this to (|68p we see that the latter is satisfied 
by the variables, (omitting arguments for brevity) 



r/ = i]--r/ + , a c = hr] + , 
f) = fj + - f)- , a c = Hfj- . 



(87) 



Inverting (|87[) we find the response substitution, 



»?+ = -£-) V+=V + 



(84) 



Unlike eqs. ([6$]). (|7D)l . eqs. ([57]). ([55]) are algebraic and 
not integral transformations. 
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. and its response transformation 



D. 



The next relation we reuse literally is eq. (|72j) . Identity 
of the linear forms in (|63[) and in (1721) requires that, 

.-, , .7 0J- -0+ ron s 
= ir?^++z a c - iip-i] - ia c , (89) 

where we applied substitution (|55|) and rearranged the 
terms. This formula hints at the variables, 



0_ -<ip + 
C = z , Vc = ip+, 

~ 0+ - 0- ~ ~ 

C = — z — , y> e = ip-, 

with the inverse substitution being, 

0+ = V'e, 0+ = "00 + K- 
i>- = 0e, 0- = 0c + ^C- 



(90) 



(91) 



Equations (|87j) . (|90|) and (|88]) . (|9Tj) are rescaled versions 
of each other. By making use of the standard chain rules 
wc find, 



5 

-l - : 

<50+ 
, <5 _ %_ 5_ 



i 5 .5 



.6 i 5 

(50o ft <^C 



, g _ , _6 i S 

50- ~~ (50c ft <5C' 



(92) 



We remind that for linear transformations with c-number 
coefficients the chain rules are identical for c-number and 
Grassmann variational derivatives, see appendix IB 41 As 
in the semirelativistic case, the purpose of additional 
phase factors in (|92[) is to emulate response substitution 
([55)1 . Indeed, the latter is found from the former by re- 
placements (|76[) ; that is yet another formula we reuse 
literally. Such replacements turn eq. (|68l) with A' - ' = 
into, 



Zc 



5?P+ ' Sip+ H- ' Sip- 



\0lpe 6( 5( 5lp e J 

proving that variables ([9U]) are indeed the right set of 
causal variables. 



E. The causal Wick theorem 

Applying transformartions (|91[) to the Hori-Wick the- 
orem (|5T|) we find the nonrelativistic counterpart of eq. 



J"Ar(0 e ,0o) = exp 



6 . 6 5 ~ S 

60c 6( S( Sip e 



X J"(0 O , 0c + K, 0c + 0e 



(94) 



The first of eqs. (1571) remains intact. It should not be for- 
gotten that eq. (|94| is subject to condition (|84|) . which 
in turn follows from (|78j) . We have thus covered all con- 
ceivable cases of interest: the real, the semirelativistic, 
and the nonrelativistic ones. 



VII. CONCLUSION AND OUTLOOK 

We have identified transformations of the functional 
(Hori's) form of Wick's theorem for operators which bring 
out its underlying physical content: propagation of free 
fields in space and time. The resulting relations are called 
causal Wick theorems. They are the main result of this 
paper. Their full power shows up in problems involv- 
ing nonlinear interactions. In particular, they will be 
instrumental in understanding electromagnetic actions 
and back-actions of distinguishable macroscopic devices, 
which will be the subject of forthcoming papers. 
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Appendix A: Formal definitions 

1. Operator orderings 

For fermionic operators, both time and normal order- 
ings are defined with the sign factor (fSU|). So, for the 
normal ordering, 



(Al) 



and so on, cf. eqs. (|47|) . Similarly, for the forward and 
reverse time orderings, 



T + 4>(x, t)$(x', t') = 6(t - t')i>(x, t)4>(x, t') 

+ s { 6(t' - t)i>(x',t')4>(x,t), 

T_i>{x, t)i>{x', t') = 9(t' - t)i>{x, t)i>{x', t') 

+ Sf6(t - t')ijj(x',t')ilj(x,t), 



(A2) 



etc. Sign rules for ordered fermionic products were dis- 
cussed at length in Ref. [3j. They reduce, in fact, to two 
statements. Firstly, under the orderings bosonic opera- 
tors behave as c- numbers (i.e., commute with any other 
operator), while fermionic operators behave as Grass- 
mann fields (i.e., anticommute with fermionic operators 
and Grassmann fields, and commute with bosonic opera- 
tors). Secondly, the order of factors "as written" acquires 
no sign factor (the rule of visual consistency) . These rules 
are evident in eqs. (|A1[) and (|A2[) . 



14 



In formulae, the T±-orderings mostly occur as double 
time ordered structures T_ ■ ■ ■ T + ■ ■ ■ . Rather than visu- 
ally keeping the operators under the T±-orderings, one 
marks the operators with the ± indices and allows them 
to commute freely. For instance, 

T^(tW(t')T + $(t")i,(t'") 

= T c 4,_(tri_(t')4, + (t")j, + {t"'), (A3) 

and so on. The ± indices serve only for ordering purposes 
and otherwise should be disregarded. Using Tc results 
in a reduction of formulae in bulk which may be truly 
dramatic (compare, e.g., expressions for response func- 
tions in J2] and in Q ) . The rule of visual consistency is 
implicit in (IA3|) . where the order of factors is identical 
on both sides. Sign rules for more complicated cases are 
evident from the example, 

TcMt)i>-(t')i>+(t") = efT(4_p)#+(t)$+V') 

= etiWT+fcttftf') = kt')T+^{t")m. (A4) 

At the last step here, the operators under the T + ordering 
were formally commuted so as to absorb the sign factor. 

The Tc-ordering is commonly visualised as an ordering 
on the so-called C-contour depicted in Fig. [TJ It travels 
firstly from t = — oo to t = +00 (forward branch) and 
then back to t — — 00 (reverse branch). One visualises 
operators marked with pluses and minuses as positioned, 
respectively, on the forward and reverse branches. The 
ordering rule for a Tp-product is consistent with the trav- 
elling order of the C-contour. 

We note in passing that, in general, equations like 
(|A2[) are mathematically troublesome. In nonlinear 
problems, Heisenberg operators are, strictly speaking, 
operator-valued singular functions. Multiplying them 
by step-functions requires specifications. Even in the 
simplest one-mode case 14j, a causal regularisation was 



necessary to fix the stochastic calculus (Ito) in phase- 
space equations. Divergences and consequent regularisa- 
tion/renormalisation procedures in quantum field theory 
[3] are of the same origin. 



2. The frequency-positive and negative parts of 
functions 



Separating the frequency-positive and negative parts 
of a function is defined as, 



f(t) = f [+) (t) + f { -Ht), 

did 



(A5) 



/ (±) (i) = J —e'^ t e(±uj)f LJl f u = j dte iut f(t). 

(A6) 

The frequency-positive and negative parts are always de- 
fined with respect to native arguments of functions, e.g., 



(±)f 



(A7) 



The W operations are conveniently expressed as integral 
transformations, 



/W(t)= / dt'S^\t-t')f(t% 



(A8) 



where 



are the frequency-positive and negative parts of the delta- 
function. Of use are the formulae, 



dtf^(t)g(t) = J dtf(t)g&\t) = J dtfW(t)g&\t), 
dtf(+\t)g( + \t) = J dtf(-\t)g(-\t)=0, 

(A10) 



where f(t) and g(t) are arbitrary functions. They may 
be verified using eq. (|A6|) and the standard equality for 
Fourier-transforms, 



For more details see Ref. [2|, appendix A. 



(All) 



Appendix B: Variational derivatives by Grassmann 
fields 



1. Definition of a Grassmann algebra 

A Grassmann algebra is a set of formal polynomials of 
the Grassmann generators 7^, k = 1, • ■ • , K. The latter 
have the properties, (fe, I = 1, • • • , K) 



Ikli = -lilk, l k = 0, 



(Bl) 



where B and F stand for arbitrary bosonic and fcrmionic 
operators. A member of the algebra is a formal polyno- 
mial, 



7 (o) 



9l ik + }^ ali'ikii 



Kk<K 



l<k<l<K 



+ &K , n'Y2-~lK- ( B2 ) 

In physical applications, coefficients 9kXi---k L i ^ = 
0, • • • , K (coefficient functions, for brevity) may be q- 
numbers as well as c-numbers. Addition of Grassmann 
polynomials is addition of coefficient functions, while 
multiplication is defined as a formal application of the 
distributive law preserving the order of factors. Consis- 
tency of formal techniques is warranted by generalised su- 
perselection rules which prohibit addition of bosonic and 
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fermionic quantities. The former (latter) is a linear com- 
bination of products of even (odd) number of fermionic 
fields and Grassmann variables. C-numbers and bosonic 
fields do not count. 

Coefficients in (|B2[) may depend on parameters such as 
space-time variables and indices (4- vector, spinor, etc.). 
So, auxiliary Grassmann variables, or Grassmann fields, 
are linear combinations of the generators where coeffi- 
cients are c-number functions Q, 



A" 



(B3) 



k=l 



Variable x has the same meaning as for the channel oper- 
ators (|47|) . Grassmann fields are fermionic. They occur 
only in the intermediate algebra and never turn up in 
observable quantitities. Properties (|B1[) extend to Grass- 
mann fields, which commute with bosonic quantities (B) 
and anticommute with fermionic ones (F): 

ip(x, t)B = Bip(x, t), ip(x, t)F = -Fip(x, t). (B4) 



2. The uniqueness theorem 

Let g(x,t) be a member of the Grassmann algebra 
where the coefficients are functions of x, t, 

K 

g(x,t) = -fez.^'*)^! "")*£• (B5) 

L=0 {k} 

To have the derivative by ip(x,t) uniquely defined we 
need the following uniqueness theorem. If 



dxdtg(x,t)Sip(x,t) = 



(B6) 



for an arbitrary Grassmann field 5tp{x,t), then g{x,t) = 
0. This theorem does not hold so far as we keep dimen- 
sionality K fixed. The obvious (and only) counterexam- 
ple is 



g(x,t) 



9X2.' 



K 



(x,t)7i • • -7A-. 



(B7) 



We therefore make K itself a variable. The limit K oo 
is not implied; K may vary but always stays finite. This 
spares us all questions of convergence. A general member 
of the algebra is then specified by (mathematically, is) a 
countable sequence of coefficient functions, by definition 
independent of K. A Grassmann field is an infinite vec- 
tor of c-number functions, also independent of K. The 
uniqueness theorem holds if "arbitrary 8ip(x, t)" includes 
arbitrary coefficient functions and arbitrary K. To show 

that the coefficient function gk^l 2 ---k L ( x ^) i s zero i it su f~ 
fices to choose, 



Sip(x,t) = Si(j L+ i(x,t)j_ 



L+l- 



(B8) 



By construction, k\ < ■ ■ ■ < < L, and the monomial 
Iki • ' 'lk L lL+\ is nonzero. Hence the coefficient at it in 
(IB6I) must be zero: 



dxdtgi L i l 2 ... kL (x,t)5^ L+1 (x,t) = 0. (B9) 

This holds with an arbitrary 5tpL + i(x, t), therefore, 

&... kL (^t) = 0. (BIO) 

3. Grassmann derivative as a variational derivative 

With the uniqueness theorem proved, derivatives by 
Grassmann fields may be introduced as variational 
derivatives, 



5F(a) = / dxdt8a(x,t) 



SF(a) 
6a(x, t) 



(Bll) 



Dependence of F(tp) on ^p{x,t) is understood as depen- 
dence of the coefficient functions of the former on those of 
the latter. This is a generalisation compared to Ref. Q 
where functionals of Grassmann variables were confined 
to formal Taylor series. Coefficient functions of 8F(ijj) 
are by definition variations of the coefficient functions of 
F(ip). The same applies to 5ip(x,t). 

If the variational derivative exists, it is unique; an 
example of a non-differentiable functional of ip(x,t) is 
ipi{x, t)j2- Straight from the definition we can prove the 
relations, 

T&tj-* Hit -'C-'W-O. <-) 

where C does not depend on ip(x) (which means that 
the coefficient functions of the former are independent 
of those of the latter). By comparing coefficient func- 
tions on both sides it is easy to show that, for any two 
diffcrcntiable functionals Fi(tp) and F^ijj), 



(B13) 



From these relations we find the laws of sum and product 
differentiation, 



8ip{x, t) 



8ip(x, t) 5ip(x, t) 



- - :F 2 [ip) + Vi^iKV) 



8tp(x,t) 8if){x,t) 



8ip(x, t) ' 
(B14) 



where rji — 1 if Fi(tp) is bosonic and r\\ = —1 if fermionic. 
Equations (|52p in the main body of the paper partly copy, 
and partly follow from, eqs. (jB12|) - (jB14|) . 



16 



4. Transformation of variational derivatives 

Assume now that we change variables in a functional, 
F(il>)=g(<p), (B15) 

where 

K(x,t,x' ,t') is an arbitrary c-number kernel. Equation 
(|B16p directly extends to variations, 

f««K(«*.*.0*K*.0. (B17) 

We assume that G{}p) is differentiable. Varying the RHS 
of (|B15|) we get, 



SQ{^p) — I dxdtSip(x,t) 



5(p(x, t) 

dxdtdx'dt'K{x, t, x', t')8ip(x', t') f . (B18) 

d(p(x, t) 

This proves differentiability of ^(tp) and the chain rule, 

£&L = [dx>dt>K(x>,t>,x,t)^L 
Sip(x,t) J 



5(p(x',t') 

^,8^x',t') sg(<p) 



dx dt 



Sijj(x,t) 5(p(x',t') 



The relation 



5ip(x',t') 
dip(x, t) 



K(x , t , x, t), 



(B19) 



(B20) 



is simply another form of (|B17|) . Under linear substitu- 
tions with c-number coefficients, variational Grassmann 
derivatives behave in a manner indistinguishable from 
conventional variational derivatives. 

It should not be forgotten that hidden behind apparent 
simplicity of eqs. (|B18p , (|B19p are infinite sequences of 
relations for coefficient functions of the Grassmann poly- 
nomials. The whole derivation hinges on the uniqueness 
theorem, which only holds with varying dimensionality 
of the underlying Grassmann algebra. 



Appendix C: Hori's form of Wick's theorem 

1. Definitions 

Let K(t,t') be a c-number kernel and ^>(t),?/>(t) a 
pair of either c-number or Grassmann fields. We wish to 
understand the action of the differential operator, 



K= [ dTdT'K(T,T')—^- 
J 01p{T 



) 5i>(r>) ' 



(CI) 



on multitime products of the fields. We also introduce 
the notation, 



*(r) = / dr K(t , i) 



5iI>(t')' 
S 

SiP(t') ' 



(C2) 
(C3) 



Variable r is thought about as time, but may in fact 
have any meaning, so far as the integration J dr and the 
symmetric delta- function 5(t,t') remain defined. This 
implies the consistency conditions, 



dr'6(T,T')X(T') = X(t), 



S 



6 



■i>(r') = 0. 



(C4) 



where X(t) may be a c-number field, Grassmann field, 
or the corresponding functional derivatives. Obviously, 
these must be defined as well, with the corresponding 
commutation/anticommutation properties. These math- 
ematical requirements are not restrictive, and the formal 
pattern of field products with contractions characteristic 
of Wick's theorem may be generalised pretty far. The 
deterring factor is ones ability to interpret the kernel 
K(t, t') as a contraction of a pair of operators, cf. eq. 
(|C12|) below. Without it, all such generalisations are 
meaningless (at least, for physics). One example when 
such generalisation makes perfect sense is the C-contour. 



2. Intermediate terms 
Using (|C4[) it is straighforward to verify the formulae, 

(C5) 



K,iI)(t)V = ^{t)KV + e f V(t)V, 
KAp{TyP = 4>(t)1CV + *(t)P, 



where V stands for an arbitrary product of the fields. The 
sign factor £f equals 1 when Tp(r), ipij) are c- numbers and 
— 1 for Grassmann variables. Now, let V n be a product 
of n factors. Applying (|C5[) recursively we see that fCV n 
is a sum of n terms, such that in each term one incidence 
of lowercase tp is replaced by uppercase 'J, 



ip(r) £f*(r), ^(r) ->• *(r) 
For example, 

K.iP{t)${t') = s^{t)${t'), 
K^{t')iI){t) = SS!(t')iP(t), 
JC$(T')<tp(T)i}(T") = *(t')V(t)^(t") 

+ e f ^(r')*(r)^(r"). 



(C6) 



(C7) 
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Terms with ^(r) or ^(r) in the rightmost position are 
zero and have been omitted. No coefficients occur, ex- 
cept that in the fermionic case terms with ^(r) acquire 
a minus. Similarly, /C 2 turns two factors uppercase, e.g., 

/C 2 ^(T)^(r')^(r")^(r'") = ^[*(t)^(t')V(t")^(t" / ) 
+ ef^(r)*(T y )V»(r")^(r"') + e^^(r)V'(T , )*(T")^(T" , )] 
= 2 [e f *(r)*(T , )V(T ,/ )^(T / ") + £ f *(r)?A(r')*(r")^(r'") 
+ ^(t)*(t')*(t")^(t"')]. (C8) 

Terms with ^(r) or ^(r) in the rightmost position are 
again dropped. In obtaining (jC8[> we also used the com- 
mutativity properties, 



/C*(r) = (r)/C, /C*(r) = *(r)/C 



(C9) 



Each term in (|CJ8[) occurs twice, depending on the order of 
replacements. Signs of terms do not depend on this order, 
so that identical terms always add up. It is clear that 
these rules will also work for KJ a applied to an arbitrary 
product of n factors. For m > n, the result is zero, 
otherwise it is a sum of terms, obtained by all possible 
sets of to replacements (|C6|) , with signs determined only 
by the latter, with the overall coefficient ml. In the below 
we call them intermediate terms. 



3. Wick terms 

Recall now that ^(t) are ^(r) are differential opera- 
tors. It is straigtforward to verify the recursion formulae, 



and 



^(t)iP{t')V = e f ^(r')*(T)P + K(t, t')V, 
%{t)iI){t')V = e f V(r')*(r)P + K (/, t)P. 



(CIO) 



(Gil) 



To visualise these relations, think about uppercase 'f's 
skipping along the product of lowercase ip's. The pairs 
4 , (t)^(t / ) and \&(t)'0(t') have no choice: they must 
(anti)commute. The pairs $?(t)x/)(t') and Vt / (r)-0(r / ) have 
a choice: they can either (anti)commute or "contract." 
Contracted pairs cancel, leaving instead c-number ker- 
nels K(t,t'). Terms where an uppercase has reached 
the rightmost position are zero, so that m uppercase 'f's 
either give rise to m contractions or "perish." (At this 
stage, we use the term "contraction" rather loosely. Its 
exact meaning is specified shortly.) Each of intermediate 
terms created by application of K. m to a product of fields 
thus results in a sum of typical Wick terms with m field 
pairs replaced by contractions. It is easy to see that any 
given Wick term is "parented" by a unique intermediate 
term (specified by left ends of contractions), and that 
all Wick terms "parented" by a particular intermediate 



term differ (distinguished by right ends of contractions). 
Hence, up to the overall coefficient of ml (and signs, see 
below), the differential operator K m applied to a product 
of fields generates a sum of all possible Wick terms with 
m contractions, each occuring only once. 



4. Sign rules 

Consider now signs of contractions. By itself, the act 
of contraction introduces no sign, cf. eqs. (|C11|) where 
the terms with c-number kernels occur with a plus. Sign 
factors due to replacements (|C6[) are absorbed by defi- 
nitions of contractions, (in Wick's original dot notation 

1) 



^•(t^V) = £f*(r)^(T') = e f K(r, t'), 
^(t')^(t) =*(r')^( T ) =K(t,t'). 
These relations are consistent with the property, 



(C12) 



(C13) 



which must hold for contractions in the true meaning of 
the term. 

The only nontrivial source of signs are therefore eqs. 
(IC10I) . (|C11I) . Sign of any Wick term may be worked 
out by formally commuting uppercase \l/'s to the left of 
lowercase tp with which they are contracted. It is obvi- 
ous that this rule coincides with the sign rule of Wick's 
theorem, which requires contracted fermionic pairs to be 
formally commuted to adjacent positions. 



5. The reordering exponent 

This way, the differential operator K m applied to an 
arbitrary product of fields generates a correct pattern of 
terms with 777 contractions with the overall factor of ml. 
In the fermionic case, anticommutativity of Grassmann 
derivatives assures correct signs of all terms. The total 
differential operator to be used in Hori's form of Wick's 
theorem is therefore found to be, 



= exp/C. 



(C14) 



Vital for consistency of this interpretation are eqs. (|C12|) . 
(|C13|) . which allowed us to redefine the arbitrary kernel 
K{t, t') as a genuine contraction of a pair of operators. 



6. Real fields 

This reasoning is easily adjusted to the real case where 
if>(r) — ip(r). We allow ip( T ) to be a Grassmann variable, 
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thus extending the Hori-Wick theorem to (hypothetical) 
Hcrmitian fermionic fields. We also assume that 

K(t,t')=£ { K(t',t). (C15) 

The K. form and the uppercase 4 r (r) are defined as, 

S 5 



K=\ j drdr' K{t,t') 
*(r) = /^,*)^. 
The rules of commutation are, 



SiP(t) Sip(T'Y 



(C16) 



K$>(t)V = V(t)1CP, 

K^(t)V = ip(r)lCV + V(t)V, (C17) 
*(r)V(r')^ = £{ip(r')^(T)V + K{t' , t)V. 

Connection to physics is established by the definition, 

( T )^' (t') = K{t', t) = E{K(t, t'). (C18) 



The only source of signs is then the last of eqs. (IC17I) . 
The rest of reasoning generalises trivially. 



Appendix D: Algebra of Green functions of a 
quantum channel 

1. The commutator and related quantities 

Up to the fact that Ar was defined with transposi- 
tion of arguments, the algebra of Green functions of a 
quantum channel is exactly that of a conventional com- 
plex field. We therefore only list formulae with terse 
comments. For more details see our papers [H, Q. For 
conventional quantised fields, many of the relations in 
this appendix appear in various textbooks on quantum 
field theory. Some ma y b e found in Schwinger's series on 
quantum field theory jlH l27j . 

All kernels we introduce are related to the commutator, 

-ihA(x, x', t-t') = (0| [4>(x, t)J(x', t')} _ et |0>, (Dl) 

where we have extended the notation with the sign factor 
to the commutator, 



[A,B] = AB - SfBA. 



(D2) 



We remind that "commutator" means commutator for 
bosons and anticommutator for fermions. Quantity (|D1|) 
is a c-number, so putting it under the averaging is a mat- 
ter of unification of all definitions. Subtracting the last 
of eqs. from the next-to-the-last one we find, 



A(x,x',t-t') = A^(x,x',t-t')+A(-\x,x',t-t'). 

(D3) 



Unlike (|D1[) , in (|59[) vacuum averaging is essential. Using 
63), 



K 

A^\x,x',t- t') = -is { ^2 



u K (x)u K (x')e 



2w K 

v K (x')v K (x)e iu ^-^ 



2lo, 



(D4) 



We see that A^ is purely frequency-positive and A^ - -* 
is purely frequency-negative. Notation Aw is consis- 
tent with these kernels being the frequency-positive and 
negative parts of the commutator. This property will be 
instrumental in expressing contractions by the the trans- 
fer functions. 



2. Response functions 



In terms of A, eqs. (|54|) and (|55[) become, 

ArO, x', t-t') = 6(t - t')A(x, x',t- t'), (D5) 

A R {x, x', t-t') = -6{t' - t)A{x, x',t-t'). (D6) 

Explicit expressions for the transfer functions are found 
using eqs. (|D4|) . 

3. Contractions 

The kernels A^) occur as contractions in Wick's the- 
orem, cf. eqs. (|55|) . The other two contractions are ex- 
pressed by them, 

A F (x, x', t-t')= 8(t - t')A(+\x, x'.t- t 1 ) 

-6(t' -t)A(-\x,x',t-t'), (D7) 

A F (x, x', t-t')= 6{t ~ t')A ( " } (a:, x, t - t') 

-6(t' -t)A^(x,x',t-t'). (D8) 
Explicit expressions again follow from (|D4|) . 



4. Response transformation of the quantum Green 
functions 



We now express contractions (|59|) by the transfer func- 
tions (O. We follow Ref. Q, the case of free complex 
fields, with two inessential amendments: there are now 
two formally independent transfer functions, and trans- 
position of arguments occuring in response transforma- 
tions in Ref. [l| is now built into the definition of Ar. 

Adding A^) to jD7| and A<+) to (|D8|) and recalling 
(|P3l). (|D% and (iDl we find, 



A F (x,x',t- t') + A l ^[x,x',t- t') = A K (x,x',t- t'), 

A F (x,x',t- t') + A {+ \x,x',t- t') = A K (x,x',t- t'). 

(D9) 
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Subtracting (|D6|) from (|D5I) yields 



and 



A R (x, x', t-t')- A R (x, x', t-t') = A(x, x',t- t'). 

(D10) 

Separating the frequency-positive and negative parts in 
this relation results in the first pair of eqs. (157)) . Substi- 
tuting them into eqs. (ID9[) and recalling (|A5[) we recover 
the second pair of eqs. (I57|) . Up to the aforementioned 
two distinctions, response transformations for the quan- 
tum channel are indeed identical to those for a complex 
field 



Appendix E: Derivation of equations (|54 [) — (|56[) 

In the first order in Hi, the Heisenberg channel opera- 
tors are given by standard expressions, 



V(x,t) =i>(x,t) - - J dt> / dx'[ip(x,t),Hi(x\t')], 

— oo 

(El) 



*(a;,t) = <ip{x,t) - - I dt' I dx'[ip(x,t),Hi(x',t')]. 



(E2) 



Consider firstly the commutator in eq. (|E1[) . 



${x,t),Hi{x',t')] 

= ijjfia' + fa'i)' - fta'tp - ctVV- (E3) 

For brevity, we write, ip = ip(x,t), ip' = ip(x',t'), and so 
on. To apply eqs. (|52p . a' in all terms must be commuted 
to the left. With eqs. (|5Tj) we have, 



ipip'a' = a'ipip' , 
tpa'ip' = Sfd'ipip' , ip' a' ip = Efc'ip'ip. 



(E4) 



In the first term, the source commutes with two opera- 
tors, causing (for fermions) two sign changes which cancel 
each other. In the other two, it commutes only once; for 
fermions, the sign change persists, resulting in the com- 
mutator turning anticommutator: 

${x,t),Hi{x',t')]=a'$,i']_ ei -a'[$,i>]_ ei . (E5) 



We use notation (ID2[) . By direct calculation, 

[i>',i>]_ ei = 0, (E6) 
and we recover eq. (|54|) . Similarly, 
$(x,t),Hi(x\t')]=<r'[lft)_ e -a'$\^], (E7) 



0. 



(E8) 



resulting in eq. (|55|) . Equations (|56l) are a consequence 
of eqs. (Eg) and (|EH) . 

Appendix F: Proof of Wick's theorem for the 
quantum channel 

We assume known Wick's theorem for the harmonic 
oscillator and for the two-level system. Namely, 



T c T(b K+ ,bt + X~,b{_) = :F N (b R ,bl):, where 
(fc) / 6 S S S 



x T( 



'ft+J 



In this relation, 



(Fl) 



(F2) 



where b K ,b\ are specified by eq. (|48l) . and b K ±(t), b K ±(t) 
are four auxiliary functional variables, c-number or 
Grassmann. The reordering form reads, 



4 fc) (/+,/+> /-./-) 



= -fe £ / dtdt'[g FK (t-t')f + (t)f + (t') 
g* FK (t' - t)f-(t)f-(t') + gW(t - ■ 



where the kernels are given by the formulae, 

gi + \t-t')=ie-^ (*-*'), 
ffFfe (i-i') = ^(i-i')e-^ (t " t ' ) , 



(F3) 



(F4) 



and the sign factor was defined by eq. ([60]) . An identical 
set of relations holds for the antiparticle operators, 



c K (t) = c K e 



ct(t) = cte*--* 



(F5) 



For fermions, the series for the exponent is truncated, 
but this does not affect the ensuing argument. 

Equations (|Fl)-(|F3) are eqs. $57$-^ downgraded to 
the nonrelativistic case for a single mode. Formulae (|F4[) 
for contractions may be verified in close similarity to eqs. 
(|59|) . The difference is that, unlike quantum channel, 
Wick's theorem for a single mode is a well established 
fact which we can rely upon. 

An arbitrary functional of the channel operators may 
be brought to the normal form by applying the Hori-Wick 
theorem (jFip modewise. The channel operators are then 
regarded not entities as in eq. (|57|) . but linear combina- 
tions of 6 K (i), c K (t), b' K (t), & K (t). The computational part 
of the Hori-Wick theorem turns into, 
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\Sb K+ ' 6b K + Sb K J 5b K . ) c \Sc K+ ' 5c K+ ' 5c K J 5c K 

I 



I b K - — b K .L ,b K ——b K ^.,c K 



, (F6) 



where 



(F7) 



Strictly speaking, eq. (|F6|) holds only if its RHS remains 
a functional of fields (|F7|) . Indeed, with eq. (IF7|) regarded 
a change of functional variables, the standard chain rules 
are, 



Sb K± {t) 
5 

5c K ±(t) 
5 

5b K ±(t) 
5 

Sc K ±(t) 



5ip±(x,t) ' 



±-Jdxu K (x) 

~ ( J < 

- — / dxv K (x) — , 

/ dxu K (x) — , 

2uj k J V 'S^ ± (x,t) 



(F8) 



h 

2uj k 



dxv K (x) 



5ip±(x,t) ' 



Differentiations in (|F6|) may therefore be rewritten in 
terms of entire fields (|F7j) . and the modewise variables 
b K (t), c K (i), b' K (t) 7 c' K (t) excluded. The key property of 
variational Grassmann derivatives is that, under sub- 
stitutions with c-number coefficients, chain rules for 
them are identical with those for conventional variational 
derivatives (appendix lB 4p . Chain rules (|F8|) thus equally 



apply to bosons and fermions. Using these relations it is 
straightforward to show that, 



E 



z. 



(fe) / _o o 6_ 

( 1 6b K + ' 5~b K+ ' 5b K ^ ' Sb K - , 
^) / 5 5 5 5 



J ° \5c K+ ' 5c K+ ' 6c K - ' 5c K 

[ 5 5 5 5 



c 



\5ifj + ' 5^ + ' 5-0- ' 5ip _ 



(F9) 



The condition, 



b K - = b K+ , b K _ = b R A 



may then be replaced by, 



Ck— ^-K— Q 



(F10) 



1>-=1>+, 1>-=1>+, (FH) 
as in (|57|) . We have thus derived the Hori-Wick theorem 
for the quantum channel from the Hori-Wick theorems 
for mode operators. 



Appendix G: Extension of the test-case formula (|63p 
to fermions 



Expanding all exponents in the intermediate formula 
in (|63[) in Taylor series we have, 



^ j^rn+m/' +n' -\-n"' ^^m+m' -\-m" -\-m ' 

ml m! ! m" \ m fN ! n ! n f ! n n ! n"' ! 

{m,n} 



xl£fT ^A F ^V7e fT ^A F ^V7efT^A (+) ^V l f^AH 6 



5ip + J V 5$- ) V &i>- 5ip + J V <*0+ Sip- , 

x (^ + V+) n (^-^-)"'(^+r?+) n "(^-r/-)""'| V;±= ^ ±=0 - (Gl) 

We use condensed notation (fH))) . Importantly, all "building blocks" in (|G1|) are bosonic: they contain an even number 
of anticommuting quantities (either zero or two). Bosonic quantities behave to a large extent as c-number ones, 

6 -a p -u (^a f ^-) = r^u^ r^A F ^ ^ (G2) 



6<lj; + J\5ip- Sip-/ \Sip- 5^-J\5^+ 5^+ 
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etc. These properties allowed eq. (IG1I) to be derived in 
the first place. 

Using (|52j) it is straightforward to derive the commu- 
tational rules, 



S S \ 6 f 5 S \ 



fjK—) (if>7]) — fjKrj + (^77) ( f/K — - ) , 



6 



(G3) 



where K{x,x',t — t') is a c-number kernel and 
ip{x, t), ip(x, t) are auxiliary variables, c-number or Grass- 
mann ones. Applying these rules to eq. (|G1|) we can 
commute all derivatives to the right; any derivative that 



reaches the rightmost position is set to zero. The result 
has the structure, 



{m} 

x (^A F?7 _)" l '[^A(+) 7?+ ] m "[r) + A(-) ?7 _] m "', (G4) 



where C mm ' m " m '" are c-number coefficients. The criti- 
cal observation is that all fermionic sign factors present 
in (|G1|) "perish" in commutations (|G3[) . Coefficients 
C m m'm"m'" therefore coincide for bosons and fermions. 
In turn, eqs. (|66[) and (IG4I) coincide up to the expansion 
of F in Taylor series. This suffices in order to extend the 
test case formula (1531 to fermions. 
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